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We study three sources of error in our calculation of Bk using HYP-smeared staggered fermions 
on the MILC asqtad lattices. These are (1) dependence on the light sea quark mass; (2) finite vol- 
ume effects; and (3) the impact of an order of magnitude increase in the number of measurements. 
Our main results are (1) the dependence on the light sea-quark mass is weaker than expected by 
naive dimensional analysis, (2) including finite volume effects in SU(2) staggered chiral perturba- 
tion theory fits leads to a very small change in Bk, of size m 0. 1 %, and (3) increasing the statistics 
on one of the coarse MILC lattices resolves a potential discrepancy with other coarse results. 
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1. Introduction 

This paper is the third in a series of four proceedings describing our calculation of Bk using 
improved staggered fermions. Here, we review some of the errors quoted in the error budgets for 
Bk given in the companion proceedings ^ and [0]. In particular, we consider the following issues: 

• The dependence of Bk on the light sea quark mass amf, 

• Finite volume effects in the SU(2) analysis of Bk; 

• The effect of increasing the number of measurements on the C4 ensemble. 

For our notations for fits, and details of the ensembles we use, see Refs. [0,1], as well as our recent 
long article [^]. 

2. Dependence of Bk on light sea-quark masses 

A year ago, we studied the light sea-quark mass dependence using five different "coarse" 
(a 0.12 fm) MILC ensembles, with results presented (using fits based on SU(2) staggered chiral 
perturbation theory [SChPT]) in Ref. [Qj. In the intervening year, we have added a second sea- 
quark mass on the fine (a w 0.09 fm) lattices (ensemble F2, with ami : am s = 0.0031 : 0.031, i.e. 
with ni( halved compared to ensemble Fl), and also increased the statistics on one of the coarse 
ensembles (C4, with ami : am s = 0.007 : 0.05). This allows us to solidify our understanding of 
the sea-quark mass dependence. This section summarizes the more extensive discussion given in 
Ref. [§. 

Both SU(2) and SU(3) analyses allow us to extrapolate our results to the physical valence d 
and s masses, assuming that the corresponding ChPT is convergent. We use the SChPT fit forms 
to correct the chiral logarithms for the fact that the sea-quark masses differ from their physical 
values. Once these logarithmic corrections have been accounted for, the remaining dependence on 
sea quark masses is analytic, and given by 

B K = b Q + b l (ami)+b 2 (am s ) + 0[(am) 2 ] . (2.1) 

Here ami and am s are the light and strange bare sea-quark masses, respectively. In SU(3) ChPT we 
have the additional relation b 2 = b\/2, while in SU(2) ChPT b 2 and b\ are unrelated. In practice, 
we can only determine b\, since, to date, all our ensembles have, for a given lattice spacing, the 
same value of am s . 

In Fig. [T], we show the dependence of Bk (after extrapolation to physical valence quark masses) 
on amf , for both coarse and fine lattices. The red and blue lines show, respectively, fits to a constant 
and a linear function. The fit parameters are given in Table |I[ We see that the slope b\ is consistent 
with zero, but has rather large errors. To compare the slopes at the two lattice spacings, we rewrite 



the expression (gjj) as 

\+b\{L P /A 2 x ) 



B k 



(2.2) 



where hp is the mass-squared of the pion composed of light valence quarks (in physical units), 
and A x is the expansion scale of ChPT, which we take to be 1 GeV. Expressed this way, the slope 
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Figure 1: Bk vs. amg for the MILC coarse ensembles (left) and for the MILC fine ensembles (right). These 
values are obtained using "4X3Y-NNLO" SU(2) fits Bk is obtained using one-loop matching and is 

evaluated at the scale jj, —2 GeV. 



coefficient b\ should be the same for both lattice spacings, up to (presumably small) discretization 
errors. Furthermore, naive dimensional analysis suggests that \b\ \ = Values for b\ are also 

given in the Table, and show that the magnitudes of the slopes are, in fact, considerably smaller 
than expected. This is not problematic, since a small value will occur some of the time. It is, 
however, serendipitous, since it reduces the uncertainty in the extrapolation to physical anif . 



a (fm) 


bo 


bx 


b\ 


0.12 


0.5624(64) 


+0.04(58) 


+0.005(74) 


0.09 


0.550(23) 


-0.93(492) 


-0.09(46) 



Table 1: Parameters of the linear fits to SU(2) results shown in Fig. [l| 

The corresponding fits for the SU(3) SChPT analysis are shown in Fig. ^| with parameters 
given in Table They give central values having somewhat larger magnitudes than the SU(2) fits, 
but they are still smaller than expected from naive dimensional analysis. Note that the slopes from 
the SU(2) and SU(3) fits need not be the same. 



a (fm) 


bo 


bi 


b\ 


0.12 


0.5691(62) 


-0.73(49) 


-0.09(6) 


0.09 


0.5448(176) 


-1.61(3.78) 


-0.15(36) 



Table 2: Parameters of the linear fits to SU(3) results shown in Fig. ||. 

Since we find no evidence for a significant dependence on m^, we assume, for our final value of 
Bk, that there is no such dependence, and do the continuum extrapolation using the mt /m s = 1/5 
lattices (including C3 and Fl). We then correct for a possible dependence on ami by including 
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Figure 2: As in Fig. |l| but for the SU(3) analysis, using N-BB1 fits (see Ref. [[|). 

in the error budgets a systematic error which is the difference in Bk between the result on the 
C3 ensemble and that obtained using linear extrapolation to amf 1 ^ D, §]. The error from an 
incorrect value of am^ hys is estimated similarly. 



3. Finite Volume Effects from SU(2) SChPT 

Finite volume (FV) dependence is predicted by ChPT. At NLO, this dependence enters through 
corrections to the chiral logarithmic functions, as follows: 



£(X)=X[log(X/^ R ) + dr(X)}, 

dl(X) _ , /v i % \ , cFV, 



£(X) 



dX 



■\og(X/^ R )-l + ^(X). 



(3.1) 
(3.2) 



Here X is the mass-squared of a pion, and /Idr is the scale introduced by dimensional regulariza- 
tion. 1 The FV corrections have the form of an image sum, 



ML ko 1*1 



8! Y (M 2 ) = 2£K (\n\ML), 

n^0 



(3.3) 



where L is the box size, and n labels the image position. For our geometries, sufficient accuracy is 
attained by keeping only spatial images. 

In our previous work, including our long article pi], we have not included these finite volume 
corrections, due to the high computational cost of implementing them. Instead, we have shown that 
the expected size of these corrections is small compared to other errors [Q]. This year, however, we 
have been able to do one-loop SU(2) SChPT fits including FV effects (keeping sufficient images 
that no approximation is made at double precision accuracy). To do so has required that we use 



The NLO result is independent of ^dr once one includes the analytic terms, and in any case /Xdr does not enter 
the FV corrections, 8fJ. 
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Figure 3: Bk{I-1 = \/a) (one-loop matched) versus Xp (mass-squared of Goldstone-taste pion composed 
of quark and antiquark of mass m x ), including 4X-NNLO fits (explained in Ref. ||). Ensembles are C3 
(left, with am y = 0.05) and SI (right, with am y = 0.018) ensembles. The red (blue) lines shows fit functions 
without (with) FV corrections included. See text for more details. 



GPUs. On a single core of the Intel i7 920 CPU (running at ~ 0.5 giga flops per core), a single FV 
fit to all ten MILC ensembles that we use takes about a week. Using an Nvidia GTX 480 GPU, by 
contrast, we have obtained a sustained performance of 67 giga flops (45% of the peak speed) [0]. 
Hence, using the GTX 480 GPU, it takes only about an hour to do the full SU(2) analysis for all 
the MILC ensembles. This is fast enough to carry out multiple fits, as is needed to estimate fitting 
errors. 

In Fig. ||, we compare the fits with and without FV corrections on a coarse (C3) and a superfine 
(S l-a pa 0.06 fm) ensemble. The red line shows the V — > °° fit, while the blue line includes FV 
corrections. The results of fits for Bk are quoted in Table [3| These correspond to the red and blue 
points in the figures, and have been obtained by setting and mt to their physical values in the 
fit function, as well as removing taste splittings and setting V — > °° in the FV form. The impact of 
using the FV fit form is very small (~ 0.1%), as can be expected from the fact that the fit functions 
hardly differ for our pion masses. It is only for smaller values of Xp that FV effects are visible. 
Note that the FV corrections have opposite signs on the two ensembles. This is due to a competition 
between two terms having different dependence on taste-splittings (and thus having different size 
on the two ensembles) [Q]. 



ID 


B K 


B K (FV) 


C3 


0.5734(46) 


0.5738(46) 


SI 


0.4914(65) 


0.4908(65) 



Table 3: Results for B^in = l/a) with and without the finite volume corrections in the fit functions. The 
fits are those shown in Fig. |[ 
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Figure 4: B K (NDR,)i = I /a) as a function of X P for the SU(2) analysis (left) and the SU(3) analysis 
(right) on the C4 ensemble. The fit types are 4X3Y-NNLO for SU(2) and D-Bl for SU(3) (see Ref. g] for 
details). Blue points and fit correspond to 1 measurement/config while red data points and fit correspond to 
10 measurements/config. 

It is well known that the one-loop prediction of FV effects gives only a semi-quantitative 
guide to their magnitude, since higher-loop effects can be important. For this reason, and to be 
conservative, we do not use the results just presented to estimate the FV systematic. Instead, we 
use the difference in Bk obtained on the C3 and C3-2 ensembles (which have different spatial 
volumes), which leads to an error estimate of 0.85% [J3j] . 

4. Effect of Higher Statistics 

Since last year, we added an additional 9 measurements on each lattice of the C4 ensemble, 
with the source timeslices chosen randomly, and different random seeds for the wall sources. In 
Ref. [f7|] we found that such additional measurements are, to good approximation, statistically in- 
dependent. This new data allows us to resolve a small puzzle we had observed in last year's data. 

In Fig. ^ we compare the low and high statistics data. The data points have shifted in a 
way which is consistent with their original statistical uncertainties, but the result is a behavior, 
particularly at low Xp, more similar to that observed on other ensembles. Examples of the behavior 
on other ensembles are shown in Refs. [l], Q|. 

This result illustrates the importance of obtaining small statistical errors, and also shows again 
how using multiple measurements on a single configuration is an efficient way of reducing errors. 
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